Continuous wavelet transform based time-scale and multi-fractal analysis of the 
nonlinear oscillations in a hollow cathode glow discharge plasma 
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Continuous wavelet transform (CWT) based time-scale and multi-fractal analyses have been car- 
ried out on the anode glow related nonlinear floating potential fluctuations in a hollow cathode 
glow discharge plasma. CWT has been used to obtain the contour and ridge plots. Scale shift 
(or inversely frequency shift) which is a typical nonlinear behaviour, has been detected from the 
undulating contours. From the ridge plots, we have identified the presence of nonlinearity and de- 
gree of chaoticity. Using the wavelet transform modulus maxima technique we have obtained the 
multi-fractal spectrum for the fluctuations at different discharge voltages and the spectrum was 
observed to become a monofractal for periodic signals. These multi-fractal spectra were also used 
to estimate different quantities like the correlation and fractal dimension, degree of multi-fractality 
and complexity parameters. These estimations have been found to be consistent with the nonlinear 
time series analysis. 



I. INTRODUCTION 

Glow discharge plasma is a typical complex medium 
exhibiting a wide variety of nonlinear phenomena such 
as chaos, noise induced resonances, self organized criti- 
cality, frequency entrainment, complex structures, etc. [TJ 
[21 HI IH 03 El [7] , which have been studied using a wide va- 
riety of analysis methods, such as the nonlinear time se- 
ries analysis, spectral analysis, etc. Nonlinear time series 
analysis like the estimation of the correlation dimension, 
fractal dimension and Lyapunov exponents reflect only 
the asymptotic and global behavior of the system but 
cannot identify the degree of chaoticity, local structures, 
etc. These estimations are also too simplistic to char- 
acterize a time series, because the estimation emphasizes 
regions of the attractor of a time series which is frequently 
visited by its orbit in the phase space [5J QUI Ell IH] • 
Further, the estimation of these exponents require very 
large data sets, and degrade rapidly with noise [8l ITB"]. 

Spectral analysis (based on Fourier techniques) has 
also been used extensively to identify the different routes 
to chaos, frequency entrainment, and many other nonlin- 
ear phenomena in plasma [TJ [51 [71 Q3]. However, the 
spectral analysis has a limitation of yielding only the 
spectral amplitude, neglecting the phase information, as 
a result of which it can not distinguish between any two 
time series of the same spectral amplitude [10 . Since the 
short time scale characteristics of a signal arc also over- 
looked, information on the time history and the local 
properties like singularity exponents and the scale shift 



(or inversely frequency shift) are lost (TTJl [TTJ [J5] . Hence 
one has to resort to a time-scale based analysis consis- 
tent with the nonlinear time series analysis, that can de- 
tect the presence of chaos, nonlinearity and other local 
structures in the system and for this purpose continuous 
wavelet transform (CWT) is one of the best candidates. 

CWT can be used to construct contours, which in turn 
are useful to study some of the nonlinear properties like 
the scale shift as well as periodic behaviour. Charac- 
terization of the weak and strong chaos is possible by 
ridge plots constructed from the CWT decomposition of 
a chaotic time series [TT]. Statistical scaling properties 
of a complex signals can be characterized by the CWT 
based multi-fractal spectra and these spectra can also 
be used to estimate several exponents like the fractal di- 
mension, correlation dimension, etc., that are useful to 
characterize nonlinear and chaotic nature of a time se- 
ries ina oina iia En Eliza E3i- 

In this paper, we have reported some results of the 
analysis of the floating potential fluctuations in a dc glow 
discharge plasma using contour and ridge plots based on 
the CWT decomposition. We have also constructed the 
multi-fractal spectra using the wavelet transform modu- 
lus maxima and from these spectra, the fractal dimen- 
sion, correlation dimension and other complexity param- 
eters have been estimated. The rest of the article is or- 
ganized as follows: we present a brief description of the 
experimental setup and results in Sections [ITJand [TlTJ re- 
spectively. Analysis of these results have been presented 
in Sections |IV[ [V] and |VI| respectively. Finally, we discuss 
and conclude the results in Section [VH] 
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II. EXPERIMENTAL SETUP 

The experiments were performed in a hollow cathode 
dc argon discharge plasma [TJ [51 [5] whose schematic dia- 
gram is shown in Fig [l] The hollow stainless steel tube 
of diameter 4 cm was the cathode and the central rod 
of diameter 2 mm was the anode. The whole assembly 
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FIG. 1: Schematic diagram of the hollow cathode discharge 
tube. Cathode is a cylindrical hollow tube and the anode is 
a rod of stainless steel. 



was mounted inside a vacuum chamber and was pumped 
down to 0.001 mbar. Then it was filled with argon gas 
up to 0.95 m bar and a discharge was struck by a dc dis- 
charge voltage. The main diagnostics was the Langmuir 
probe used to monitor the floating potential fluctuations. 
The typical plasma density was about 10 7 cm -3 and the 
electron temperature estimated was about 3 — 4 eV. 



III. FLUCTUATION RESULTS 

The formation of the plasma depends on two main pa- 
rameters: the discharge voltage (DV) and the neutral 
gas pressure. An anode glow was formed around the 
anode when a discharge was struck above a threshold 
neutral pressure and DV. It was observed that at the 
initial stage of the discharge, the anode glow was large 
in size and decreased with increase in DV, until it got 
extinguished beyond a certain DV pQ. At the instant 
of formation of the anode glow, the associated fluctua- 
tions in the floating potential were chaotic which became 
quasi-periodic with increase in the DV and finally van- 
ished to a fixed point when the anode glow disappeared. 
At this time the plasma inside the hollow cathode be- 
came brighter and filled up the entire plasma column. 
Details of the anode glow and the associated instabilities 
have been discussed in Ref [T]. For the present paper, 
the floating potential fluctuations were obtained keeping 
pressure constant at 0.95 m bar. At this pressure the 
discharge was initiated at ~ 283 V and simultaneously 
irregular oscillations in the floating potential were ob- 
served [Fig^a) (left panel)]. With increasing DV, the 
regularity of the oscillations increases and final form of 
the fluctuations just before their disappearance was the 
relaxation type of oscillation as shown in Fig |2^h—i) [left 
panel] . 



IV. CWT ANALYSIS OF THE FLUCTUATIONS 

The continuous wavelet transform (CWT) is used to 
decompose a signal using wavelets, i.e., into small oscil- 
lations that are highly localized in time and for a signal 
4>(t) it is defined as [22] : 



Wy(a,T) = / <t>(t)V a (t-T)dt 



(1) 



where, <j)(t) is the signal and ^f a (t) is an oscillating func- 
tion that decays rapidly with time and is termed as 
wavelet, a and r are the scale and temporal propaga- 
tion parameters respectively. Scale 'a' can be consid- 
ered to be the inverse of frequency and so one can es- 
timate the frequency (F a ) corresponding to a scale 'a' 
for a particular wavelet using the relation F a = -^t (24] . 
where F c is the center frequency of the analyzing wavelet 
and A is the sampling period of the signal. So a signal 
can be decomposed in time-scale plane or time-frequency 
plane using above scale-frequency relation. In this pa- 
per, we have used the time-scale decomposition using the 
Daubechies orthogonal wavelet to detect the singularities 
present in the time series as it has a slight asymmetric 
structure, uses few coefficients and is a good represen- 
tation of low-order polynomials [25] 126] . Fig [2ja'— i') 
shows the contours of W^(a, r) for the signals shown in 
Figs [2ja — i), in the time-scale plane. Though one can 
see that the oscillations go through different types of 
behaviour with increasing DV, changes are much more 
clearer in the CWT contour plots shown in Fig [2] (right 
panel). Undulating nature of the contours in Figl^Ta'— b') 
and Fig[2jf' g') is due to the scale shift with time which 
is one of the signatures of the presence of nonlinearity 
in the system [27] . Periodic behaviour of the fluctuating 
oscillations of Fig[2^c— e) is also prominent from the pe- 
riodic contours [Fig [2jc' — e')]. Though the contour plots 
of the relaxation oscillations [Fig [2 (h— i)] have no such 
distinct structures, still one can find differences from the 
other two cases. Importance of the CWT based contour 
plots is that they can be used to identify the nonlinear 
nature of a system, and this identification is very essen- 
tial for the nonlinear time series analysis. The CWT 
coefficients (VF$(a,r)) are also useful to devise another 
plot termed as a ridge plot [TT] which can be exploited 
to analyse chaotic data, and this has been presented in 
the next section. 



V. RIDGE ANALYSIS 

The line joining the maximum of the CWT coefficients 
along the scales 'a' is called a ridge plot which can be used 
to detect the presence of chaos and degree of chaoticity 
or periodicity in a system [11] . Here, we have consid- 
ered only the prominent ridges to construct the plots. In 
Fig [3](a) the presence of the broken ridges at the higher 
scale of 200 and also at the lower scale of about 30 indi- 
cates that at the initial stage of DV (283 V) the system is 
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FIG. 2: (Left panel) Floating potential fluctuations (a— i) and (right panel) their corresponding contour plots of continuous 
wavelet transforms (a'-i') are shown at 0.95 m bar for different voltages: (a) 283 V; (b) 284 V; (c) 286 V; (d) 288 V; (e) 289 
V; (f) 290 V; (g) 291 V; (h) 292 V; (i) 293 V. Nonlinear scale shifts are clearly seen in a' and b'. c'— e' is almost periodic and 
with further increase in DV the signals exhibit nonlinear scale shift before becoming periodic relaxation oscillations (V— i'). 



very chaotic [IT]- At 284 V [Fig^b)] the system is still 
chaotic but the discontinuous ridges are seen only at the 
lower scales and chaos is not so strong and these results 
are consistent with the nonlinear time series analysis [1] . 
With increase in the DV, the fluctuations [Fig[3jc— e)] be- 
came almost regular and their ridges are almost constant 
with time. The ridge plots of the fluctuations shown in 
Figs|2jf) and (g) also show broken ridges indicating reap- 
pearance of the chaos. Appearance of a chaotic window 
after the periodic case is always possible and is consis- 
tent with previous analysis presented in Ref pQ. Fig 3^h) 
shows discontinuous ridges but of slightly longer dura- 
tions and finally Fig [3](i) shows a continuous ridge in- 
dicating periodic signals with breaks at the fall phase 
of the relaxation oscillations at the lower scale of about 
300. Here it also shifts between two scales one at the 
top and the other at the bottom indicating the nonlinear 
frequency shift. 



VI. MULTIFRACTAL SPECTRUM 

In order to understand the scaling properties of the 
signals we have constructed multi-fractal spectrum us- 
ing the wavelet coefficients of the CWT decomposition. 
The multi-fractal spectrum (D(a)) of a signal can be es- 
timated using the wavelet transform modulus maxima 
method(WTMM) [IB [HI [22 I2H] as follows: 

The partition function Z(a, q) of moment q for a partic- 



ular scale 'a' is defined as a sum of the wavelet transform 
modulus maxima at that particular scale is given as 



Z(a,q) 



E 

{ti(a)}i 



|W«(o,t i (o))|«~o T W 



(2) 



where, W^(a,U(a)) is the wavelet modulus maxima and 
q is the moment considered. From r(q) we can estimate 
the singularity spectrum using following relations |22j 



(3a) 



D(a) 



drjq) 1 
dq 2' 
min(g((a+l/2)-T(g)), (3b) 



We have estimated the multi-fractal spectra for all the 
nine signals shown in Fig[2^1eft panel) at different DV. In 
Fig [4] black dots (— • — ) show the typical multi-fractal 
spectra (D(a) vs a) at 283V for the initial stage of dis- 
charge and uptriangle (—A—) for the fluctuations at 292 
V just before the stable state. These plots show that the 
signals at the initial stage of discharge is multi-fractal 
in nature and tends to become a mono-fractal with in- 
crease in DV (for quasi-periodic signal) as the width and 
height of the spectrum have decreased considerably. For 
the periodic signals [Figs [2|d, e and i)], we obtained a 
mono-fractal behaviour. 

The correlation dimension (D corr ) has been estimated 
from their multi-fractal spectra using the relation D q= 2 — 
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FIG. 3: Ridge plots along the scales of the contour plots shown in Fig [2] (right panel). 



1.0 - 



0.8 - 



0.2 - 



0.0 




1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 

a 

FIG. 4: — • — shows the multi-fractal spectra at the initial 
stage of discharge (283 V). —A— at 292 V, i.e., just before 
the steady fixed point is reached and mono fractal in nature. 



2a g= 2~ f{ot q= 2) [13], has been shown in Figp^a) by black 
dots (—•—). It is seen that the D corr is nigh for the 
initial stage of discharge and decreases with DV except 
for 290 V and 293 V. The increase in D corr at 290 V is 
probably due to the appearance of a chaotic feature at 
an intermediate stage of the periodic signal. 

For comparison, D corr using Grassberger-Procaccia al- 
gorithm pQ, has been shown in the same plot (open cir- 
cle) and both exhibit a good agreement. There are slight 
discrepancies observed at 283 V, 286 V and 293 V. The 



discrepancies at the first two values can be possibility 
a result of the presence of modes at significantly differ- 
ent scales. This would mean that since wavelet trans- 
form is based on the multi-resolution analysis of vary- 
ing frequency and time resolutions at differing scales, 
the WTMM estimates joining maxima at different scales 
might incur errors in the estimated values of D corr . How- 
ever, the general trend of decrease in D corr , from the 
two techniques, is quite similar in nature. The increase 
in D corr at 293 V may be due to nonlinear nature of 
the relaxation oscillations. The degree of multi-fractality 
(/?) defined as the difference between the maximum and 
minimum values of a and which also statistically gives 
the range of scale invariance |29j is shown in Fig^c). 
Since (3 is also a measure of complexity, it is seen that 
the complexity [Fig [5^ a)] decreases with increase in DV. 
Fractal dimensions, i.e., f(a q= o) of the singularity sup- 
port of the signals at different the DV have been shown 
in Fig |5jc) and they decrease with increase in DV. Fi- 
nally we have estimated another complexity parameter 
(P c ) which is high for the complex signal [30], and is 
defined as P c = jj^Soyj H FW H M , where, h max is the 
singularity exponent at the fractal dimension (f(a q= o)) 
and UpwHM is the full width at half maxima of the spec- 
trum. Fig[5jd) shows the behavior of P c at different DV. 
At the initial stage of the discharge P c is high and de- 
creases at 288 V and 289 V, which corroborates well with 
estimated D corr [Fig[5|a)]. P c increases at 290 V where 
the signal is complex [Fig (5^c)], and for the last case we 
have high P c which may be due to the nonlinear effect of 
the relaxation oscillation [Fig^i)]. 
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FIG. 5: (a) Correlation dimension (D corr ) obtained from the multi-fractal spectra (— • — ) and D corr from Grassberger-Procaccia 
algorithm (open circle); (b) degree of multi-fractality (/3); (c) fractal dimension (Df)', and (d) complexity parameter (P c ) for 
different DV. 



VII. DISCUSSION AND CONCLUSION 

Glow discharges are simple systems but exhibit exotic 
features depending on the configuration, discharge pa- 
rameters, etc. The complexities in the plasma dynamics 
arise from many degrees of freedom like different sources 
of free energy, various types of wave particle interaction 
and many other instabilities. As the frequencies of the 
instabilities presented in this paper, were within the ion 
plasma frequency and the presence of anode glow and 
relaxation oscillations have been attributed to the pres- 
ence of double layers, these instabilities could be related 
to ion acoustic waves in the presence of ionization insta- 
bilities and also anode glow related double layer [T]. Such 
double layer associated ionization instabilities have also 
been observed before [31]. Different analysis tools have 
proved to be very useful in exploring the various possible 
linear and nonlinear features of these instabilities. 

Here, we have carried out wavelet based time-scale and 
multi-fractal analysis of the floating potential fluctua- 
tions of the glow discharge plasma for extracting some 
of the nonlinear features. From the contour plots of 
the CWT coefficients, we have detected scale shift (or 
inveresly frequency shift) with time in the glow dis- 
charge plasma. Though theoretical prediction of such 
nonlinear phenomenon has been reported in plasmas 
[32] 133] 134] 135] . there are very few experimental ev- 
idences on this [3B], and we have identified this phe- 
nomenon using wavelet techniques. The scale shifts can 
occur if the coherent modes like the ion acoustic wave or 
a double layer propagates through an ion acoustic turbu- 



lence [3U [35] • The ridge patterns of the plasma signals 
based on CWT decomposition reveal the chaotic features 
of the system and these finding agree well with our earlier 
analysis using the nonlinear time series analysis pQ. The 
estimated correlation dimension, fractal dimension, and 
complexity parameter from multifractal spectra are also 
in close agreement with the nonlinear analysis [1.. With 
increase in the DV, the multi-fractal behaviour changes 
to a monofractal as the irregular oscillations become pe- 
riodic with increase in the DV. 

Finally, we feel, the analysis of chaotic data using the 
contours, ridge plots and multi-fractal spectra, based on 
the CWT decomposition, are very useful because these 
do not require a long data length. However one should 
be cautious about the choice of the appropriate wavelet. 
Whereas, the estimation of the correlation dimension and 
Lyapunov exponents estimated from the nonlinear time 
series analysis require a long data length, and these esti- 
mations depend also on the proper choice of the embed- 
ding dimension, time lag, etc. 



acknowledgement 

We acknowledge the constant support and encourage- 
ment from the Director, SINP. We would also like to 
thank our colleagues of the plasma physics division for 
their help during the experiments. One of the authors 
(MN) would like to acknowledge the support of Amit 
Apte at TIFR. 



G 



[1] Md. Nurujjaman, Ramesh Narayanan, and A. N. Sekar 
Iyengar, Chaos 17, 043121 (2007) and references there 
in. 

[2] Md. Nurujjaman, A.N. Sekar Iyengar, and P. Par- 
mananda, Phys. Rev. E 78, 026406 (2008). 

[3] Md. Nurujjaman, and A.N. Sekar Iyengar, Phys Letts A 
360, 717 (2007). 

[4] Lin I and Jeng-Mei Liu, Phys. Rev. Lett. 74, 3161 (1995). 

[5] A. Dinklage, C. Wilke and T. Klinger, Phys. Plasmas 6, 
2968 (1999). 

[6] T. Klinger, F. Greiner, A. Rohde, A. Piel, and M.E. 

Koepke, Phys. Rev. E 52, 4316 (1995). 
[7] R.H. Abrams. E.J. Yadlowsky, and H. Lashinsky, Phys. 

Rev. Lett. 22, 275 (1969). 
[8] Julien Clinton Sprott, Chaos and Time-Series Analysis 

(Oxford University Press, Oxford, 2003) p. 123. 
[9] Holger Kantz and Thomas Schreiber, Nonlinear time se- 
ries analysis (Cambridge University Press, Cambridge, 

2004) Chapter 5-6. 
[10] Junfeng Sun, Yi Zhao, Tomomichi Nakamura, and 

Michael Small, Phys. Rev. E 76, 016220 (2007). 
[11] C. Chandre, S. Wiggins, and T. Uzer, Physica D 181, 

171 (2003). 

[12] Sprott, J. C. and Rowlands, G., International Journal of 

Bifurcation and Chaos 11, 1861 (2001). 
[13] Gabriel B. Mindlin and R. Gilmore, Physica D 58, 229 

(1992). 

[14] R. W. Boswell, Plasma Physics and Controlled Fusion 

27, 405 (1985). 
[15] T. Higuchi, Physica D 31, 277 (1988). 
[16] A. Arncodo, Y. d' Auubenton-Carafa, E. Bacry, P.V. 

Graves, J.F. Muzy, and C. Thermes, Physica D 96, 291 

(1996). 

[17] Alexander Silchnko and Chin-Kun Hu, Phys. Rev. E 63, 
041105 (2001). 

[18] Ashvin B. Chhabra, Charles Meneveau, Roderick V. 
Jensen, and K. R. Sreenivasan, Phys. Rev. A 40, 5284 
(1989). 

[19] Zbigniew R. Struzik and Arno P.J.M. Siebes, Physica A 

309, 388 (2002). 
[20] J.F. Muzy, E. Bacry, and A. Arneodo, Phys. Rev. Lett. 

67, 3515 (1991). 
[21] J.F. Muzy, E. Bacry, and A. Arneodo, Physical Review 

E 47, 875 (1993). 
[22] Stephane Mallat, A Wavelet Tour of Signal Processing 

(Academic Press, Burlington, 1999) p. 79. 
[23] H. Grussbach and M. Schreiber, Phys. Rev. B 51, 663 

(1995). 

[24] scal2frq is a standard MATLAB function to compute 
frequency from scale. 

[25] Odim Mendes Jr, Margarete Oliveira Dominguesb, Aracy 
Mendes da Costa, and Alicia L. Cilia de Gonzalez, Jour- 
nal of Atmospheric and Solar- Terrestrial Physics 67, 
18271836 (2005). 

[26] Margarete Oliveira Domingues, Odim Mendes Jr., and 
Aracy Mendes da Costa, Advances in Space Research 
35, 831842 (2005). 

[27] Rick Lind, Kyle Snyder, and Marty Brenner, Mechanical 
Systems and Signal Processing 15, 337 (2001). 

[28] Plamen Ch. Ivanov, Luis A. Nunes Amaral, Ary L. Gold- 
berger, Shlomo Havlin, Michael G. Rosenblum, H. Eu- 



gene Stanley, and Zbigniew R. Struzik, Chaos 11, 641 
(2001). 

[29] Xia Sun, Huiping Chen, Yongzhuang Yuan, and Ziqin 

Wu, Physica A 301, 473 (2001). 
[30] Yu Shimizu, Markus Barth, Christian Windischberger, 

Ewald Moser, and Stefan Thurnerb, Neurolmagc 22, 

1195 (2004). 

[31] James C Johnson, Nicola D'Angelo, and Robert L Mer- 

lino, J.Phys.D:Appl.Phys 23, 682 (1990). 
[32] John M. Dawson, Phys. Rev. 113, 383 (1959). 
[33] G. J. Morales and T. M. O'Neil, Phys. Rev. Lett. 28, 417 

(1972). 

[34] J.Vaclavik and K.Appert, Phys.Fluids 23, 1801 (1980). 
[35] O.Ishihara and A.Hirose, Phys.Fluids 10, 2429 (1984). 
[36] M. E. Koepke, A. Dinklage, T. Klinger, and C. Wilke, 
Phys. Plasmas 8, 1432 (2001). 



